Atomistic modeling of the electrostatic and transport properties of a simplified 

nanoscale field effect transistor 
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A first-principle model is proposed to study the electrostatic properties of a double-gated silicon 
slab of nano scale in the framework of density functional theory. The applied gate voltage is 
approximated as a variation of the electrostatic potential on the boundary of the supercell enclosing 
the system. With the electron density estimated by the real space Green's functions, efficient 
multigrid and fast Fourier Poisson solvers are employed to calculate the electrostatic potential 
from the charge density. In the representation of localized SIESTA linear combination of atomic 
orbitals, the Kohn-Sham equation is established and solved self-consistently for the wavefunction of 
the system in the local density approximation. The transmission for ballistic transport across the 
atomic silicon slab at small bias is calculated. The charge distribution and electrostatic potential 
profile in the silicon slab versus the gate voltage are then analyzed with the help of the equivalent 
capacitive model. Quantum confinement and short gate effects are observed and discussed. 

PACS numbers: 



I. INTRODUCTION 

The size of Si based metaLinsulator-semiconductor 
(MIS) devices has been shrunk aggressively to nanoscale 
nowadays [ll] in design of very large scale integrated cir- 
cuits. The quantum-mechanical nature of transport be- 
comes inevitably prominent and there are extensive ef- 
forts devoted to modeling quantum behaviors in nano- 
MIS devices based on the effective mass approximation 
or the empirical tight-binding models. M 0, [E H, S 
H, Besides the traditional top-down technologies, 

bottom-up techniques has also been developed to fabri- 
cate Si nanostructures such as nanowires for future nan- 
odevices. [ll| As a result, the first-principle study on 
electronic and transport properties of atomic Si systems 
becomes attractive. [H, [3| 

Due to the extensive application of Si based field ef- 
fect transistors, some efforts have been carried out to 
the first-principle understanding of their behaviors in the 
past years. Evans et al. [13] have studied the effect of 
Si-Si02 interface roughness on the electron mobility in a 
Si based MIS structure by describing the roughness by 
first principles. Based on the above method, Hadjisav- 
vas et al. ,15| have proposed an explanation for electron 
mobility enhancement in strained MIS field effect transis- 
tors (FETs). Fonseca et al. _16J and Liu et al. have 
established a two-terminal system of Si-Hf02-Si and Al- 
Si02-nSi respectively described by a self-consistent den- 
sity functional theory (DFT) and studied the transport 
properties through the Hf02 and Si02 slab in the frame- 
work of the noncquilibrium Green's function (NEGF). In 



this way, they have estimated the leakage current through 
the ultra-thin oxide barriers of MISFETs as a function 
of the voltage drop. Landman et al. [l^ have stud- 
ied the quantum transport through short Si nanowires 
passivated by hydrogen atoms and attached between Al 
electrodes by solving the eigenchannels of the scatter- 
ing states in the framework of the density functional 
method with plan wave basis. Similarly, Ng et al. [Toj 
studied the quantum transport through hydrogenated Si 
nanowires with Li electrodes employing the NEGF-DFT 
approach integrated in the Atomistix ToolKit (ATK) 
package. Fernandez-Serra et al. [l^] have investigated 
the effect of doping on electronic quantum transport in 
Si nanowires by assuming the same material (Si) for elec- 
trodes and the device regions. Later on, Markussen et 
al. [2]| have continued the work for longer wires and 
studied the crossover from ballistic to diffusive transport 
based on the DFT and a recursive Green's function ap- 
proach. Besides the above two-terminal models, three 
terminal models have also been used for Si systems to 
take into account the gate effect. Using the NEGF-DFT 
approach with the help of ATK package, Dai et al. [2^ 
have investigated the transport properties of a Si4 cluster 
sandwiched between two Al electrodes. By assuming a 
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constant potential shift in the molecular region [2 
when a gate voltage is applied to the third terminal, they 
have observed charge transfer from the molecular and 
transconductance oscillation. 
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The aim of this work is to go further from small 
three-terminal molecular systems to more realistic multi- 
terminal MIS devices of nanoscale using the NEGF-DFT 
approach. However, a rigorous treatment [25] of realistic 
multiterminal systems is computationally too expensive. 
Fortunately, particle exchange between the gate and the 
channel is usually very limited as they are separated by 
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the insulator slab and the gate terminal works as a con- 
troller of electrostatic potential in the device region. To 
assure the simulation feasible with the present computa- 
tion capability while catching the most important char- 
acteristics of nano MISFETs, we use a simplified model 
and focus ourselves on the gate effects on the electrostatic 
and transport properties under a small source-drain bias 
(in the linear region). As we know, in a real MIS device, 
how the gate voltage is distributed into the channel of 
the device or the electrostatic property of the device as 
a MIS capacitor is critical in determining the character- 
istics of the whole device. [1^ For example, double-gate 
(DG) structure has been considered as the most promis- 
ing device geometrical structure for nano MIS devices 
due to its good electrostatic integrity [l^, [IHl- In our 
model, we will take the effect of gate voltage distribu- 
tion into account instead of simply adding a constant 
shift to the electrostatic potential. The gate voltage that 
controls the boundary condition and the corresponding 
variation of potential inside the device is self-consistently 
determined by Poisson equation. 

The rest of the paper is organized as follows. In Sec. 
II, we define the device model and briefly describe the 
NEGF-DFT formalism and the computation method. In 
Sec. Ill, the result for electrostatic properties, charge 
and potential, are reported and analyzed with the help 
of density of states (DOS) and transmission spectra. Fi- 
nally, a brief summary is given in Sec. IV. 



II. MODEL AND METHOD 
A. Device Model 

We consider a prototype DG MISFET as schematically 
shown in Fig. [T](a). The device is made from a Si slab 
with thickness of nm order in the x direction. The chan- 
nel (z) direction is chosen along Si crystal [001] orienta- 
tion and the interfaces along (100) which is the preferred 
direction due to lower interface defect density and higher 
mobility. H atoms are used to passivate the dangling Si 
bonds on the surfaces. The H-Si bond length is relaxed 
and has a value of 1.49 A while the Si-Si bond length 
is 2.35 A. Further on either side of the Si slab a metal 
gate of variable length along z direction is located a few 
angstroms away from the H atom sheet, as indicated by 
the thick solid bars in Fig. [BJa). Note that here we in- 
sert the vacuum slabs, instead of oxide slabs, between the 
gates and the Si slab as the insulator. This is one of the 
simplifications introduced in this model in order to re- 
duce the computation load. Since the dielectric constant 
of vacuum («;„ = 1) is much lower than that of oxide (e.g. 
Hox — 3.9 for Si02), we expect a lower insulator capac- 
itance (between the gate and the channel) in our model 
and, hence, a larger gate threshold than that for devices 
with oxide insulator slabs. 

Four different model structures of total length 21.72 A 
are studied in this paper and they are denoted as Smn 



for models with short gate (length 5.43 A) and Lmn for 
the model with long gate (in the full central region of 
length 10.86 A). The index m (n) is used to distinguish 
models with different thickness of the Si slab (vacuum 
slab). The Si slab in the models is m unit-cells thick 
(5.43 and 10.86 A for m = 1 and 2 respectively). The 
corresponding vacuum slab thickness is 1.9 and 2.9 A for 
rt = 1 and 2 respectively. 

Another simplification is to approximate the gate effect 
with a boundary condition based on the following phys- 
ical consideration. When we apply a small voltage drop 
between the Si slab and the metal gates we shift the Fermi 
energy in the latter by letting the former grounded. Due 
to the high electron density in metal, the electrostatic 
potential in the gates will shift in parallel with the gate 
Fermi energy or the voltage drop. As a result, the role of 
the metal gates on the channel is to control the electro- 
static boundary of the supercell as enclosed by the solid 
lines in Fig. [TJa). Understanding this gating mechanism, 
we can exclude the gate materials from our model system 
(inside the supercell) while taking its most important ef- 
fect into account by shifting the boundary condition as 
further explained in next section. 

A third simplification is to use the same material for 
the source/drain electrodes and the channel. In real MIS- 
FET devices, the drain and source electrodes are doped 
and the contact between the electrodes and the channel 
adds to the complexity of the system and makes the re- 
sult difficult to analyze. In this work, we focus on the 
gating effect in the channel region and it would be more 
practical to handle a simplified system. This simplifica- 
tion is expected not to affect the conclusion when the 
bias is low and has been used by other authors. [lO] 

After the above simplifications, in our model, the 
double-gated MISFET is a Si slab with known electro- 
static boundary and can be studied by the NEGF-DFT 
method. The system is composed of a hydrogen passi- 
vated Si slab and two vacuum slabs. Along z axis, as 
separated by the vertical dash lines in Fig. [TJa) , the sys- 
tem is divided into three regions: the left electrode (L) 
the central region, or the channel region (C), and the 
right electrode (R). In z direction, the electrostatic po- 
tential of a Si slab without gate is used as the boundary 
condition of the electrostatic potential on the left and 
right boundaries of the supercell. This implies that the 
left and right electrodes are assumed semi-infinite. In y 
direction, the periodic boundary condition is applied to 
assume an extending system in this direction. In x direc- 
tion, the boundary is variable depending on the system 
size and the gate voltage and will be specified in the cor- 
responding context. Furthermore, we assume a uniform 
Fermi level inside the supercell. This approximation is 
justified in case of small bias between the source and 
drain electrodes because the gate voltage induced varia- 
tion of Fermi energy happens mainly in the vacuum slab 
and the density of states there is negligible. 
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FIG. 1: (a) Atomistic scheme of a model double-gate MIS- 
FET enclosed in one supercell of length 21.72 A (four Si unit 
cells). The system has the left electrode (L), the central re- 
gion (C), and the right electrode (R) separated by the vertical 
dashed lines at z = 5.43 and 16.29 A. The atomic positions 
projected to the x-z plane are illustrated by the big (Si) and 
small (H) spheres. The short thick bars (marked by the gate 
voltage Vg) indicate the position of the gates. The gate re- 
gion is the supercell region between the two bars with the gate 
length along z. (b) The two equivalent capacitive circuits for 
the double gate MISFET. The surface potential Vs is the av- 
eraged potential change on the outer surfaces of the H slabs 
(indicated by the horizontal dashed lines in (a)) in the gate 
region when a gate voltage is applied. The central potential 
Vc is the averaged potential change on the central y-z plane 
in the gate region. 



density p{r). The Kohn-Sham (KS) equation [30| is then 
established for the atomic system in the representation 
{ipi} with the corresponding KS Hamiltonian operator 
expressed in matrix form as H[p{r)]. In extending sys- 
tems, we use the supercell technique with the periodic 
boundary condition to define and solve the KS equation 
in the whole real space. The dimension of the supercell is 
chosen larger than the screening length of the KS poten- 
tial in any direction so the real system is well mimicked by 
that inside the supercell. In open systems, the wavefunc- 
tion differs from that in closed system and this is taken 
care of by using the Green's function of a semi-infinite 
lead ^38*1 as the boundary condition for the Green's func- 
tion inside the supercell. With the KS Hamiltonian H we 
calculate the retarded (advanced) Green's function (GF) 
G''^(^)(e), Keldysh GF G<(e) and retarded (advanced) 
self-energy (e). The density matrix p{r) then reads 

MMm 



p 



'2^ 



j deG<{e). 



(1) 



The final density matrix is self-consistently reached from 
an initial guess, for which we use the neutral atom den- 
sity matrix obtained by assuming no interaction between 
atoms. In the linear response region where the bias 
between the electrodes is small, the Fisher-Lee relation 
pot connects the Green's functions with the transmission 
coefficients. With the density matrix and the Green's 
functions known, the transmission T(e) for the quantum 
transport between the electrodes is then given by [25j 



r(e) = Tr[riG«rflG^], 



(2) 



where F^, F^ are line- width functions of left and right 
electrodes which indicate the corresponding coupling 
strength between the electrodes and central region. The 
linear zero-temperature conductance is determined by 
the transmission at Fermi energy 
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G = {2e^/h)T{EF). 



(3) 



To analyze the electrostatic and transport properties 
of the atomic scale MISFET, we utilize density func- 
tional theory (DFT), [2^ |33] which converts the many- 
body system into a single-particle system, and the two- 
terminal nonequilibrium Green's function (NEGF) tech- 
nique [2^, (3l| , which takes care of the coherent transport 
through the system between the source and the drain 
electrodes in Landauer picture [H, [H, [H, [3^. [ssl . With 
standard norm-conserving pseudopotentials [36| to de- 
scribe the effective interaction of the valence and core 
electrons, we construct a localized SIESTA linear combi- 
nation of atomic orbitals (LCAO) basis set {i/'i} [IJ] as 
the representation for expanding the electronic wavefunc- 
tions. According to the Hohenberg-Kohn theorem, [2^ in 
the Born-Oppenheimer approximation, the fully interact- 
ing electron problem for the ground state can be mapped 
into a variational problem in terms of a single-particle 



C. Computational Method 

Atomistix ToolKit (ATK) [41|] for two-probe systems 
together with a multigrid Poisson solver [12, is used 
to carry out the numerical calculation. We at first apply 
the periodic boundary condition in the x and y directions 
and solve the corresponding electron density po{x,y,z) 
and electrostatic potential Vb(a;, y, z). VQ{x,y,z) on the 
boundary is denoted as the boundary condition at Vg = 
0. The boundary condition at finite Vg is obtained by 
shifting Vo{x, y, z) with Vg on the gate but keep the same 
as Vo(a;, y, z) on the electrodes. On each boundary region 
between a electrode and a gate the potential shift de- 
creases linearly from Vg to zero. The modified boundary 
condition is then applied to the system to solve p(x, y, z) 
and V(X,y, z) at arbitrary Vg. In the calculation we use 
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the SZP, single C valence s and p orbitals (s orbitals) plus 
single ^-polarization d (p) orbitals for Si (H), real space 
SIESTA LCAO basis set [H, M, EI and the standard 
nonlocal norm- conserving pseudopotential [36l | from the 
database in ATK. The reason why we use the SZP basis 
set instead of a simpler one is because the result presented 
in this paper converges for SZP and more complete basis 
sets. The local density approximation (LDA) with the 
Perdew-Zunger parametrization [4^ of the correlation 
energy for a non spin-polarized homogeneous electron gas 
[4^ is used for the exchange-correlation functional. The 
mesh cutoff is set to 4348.48 eV corresponding to a grid 
size of 0.092Ax0.092Ax0.092A. The MonckHorst-Pack 
k-point grid is (1, 1, 100) along {x,y, z) directions. 



III. RESULTS AND DISCUSSIONS 

To express a physical property A in our microscopic 
model on the same footing as its counterpart in tradi- 
tional macroscopic model, we measure it from its value 
at = and use its average over the y-z range inside 
the central unit ceU (8.15 A < z < 13.58 A), (^)y,z, 
in the following analysis. The total net charge Q is the 
total charge transferred into the supercell when Vg is ap- 
plied. Note that there may be a shift between Vg in 
our model and the observed in experiments because i) 
the real boundary condition at = may differ from 
the periodic (bulk) boundary condition used in the cal- 
culation, ii) the boundary potential shift when a finite 
Vg being applied is determined selfconsistently and can 
vary in different environments, and iii) no workfunction 
is specified for the gate metal. 

Equivalent capacitive circuits are widely used in quali- 
tative analysis of electrostatic properties of MISFETs in 
the literature. Our system is macroscopically symmetric 
about the y-z (top-down) and x-y (left-right) planes and 
three capacitors, Ci, Csz and Csx may be used for this 
purpose. As shown in Fig. [Hb), left circuit, Ci accounts 
qualitatively for the average voltage drop across the in- 
sulator slab while Csx for the one across the Si slab in x 
direction and Csz in z direction. For long channel as usu- 
ally exists in traditional DG MISFETs, one-dimensional 
approximation is valid and only Csx is important for the 
electrostatic analysis. For short devices, however, the 
boundary at the ends of the channel also plays an impor- 
tant role in determining the potential profile in the cen- 
ter of the channel and the effect of Csz should be taken 
into account. Because only the linear region is concerned 
in this paper, we can also use a single capacitor Cs to 
describe the electrostatic properties of the Si slab with 
Cs = CsxCszl [Csx + Csz) ^ Csx, as illustrated in the 
right circuit of Fig. [ijb). The total net charge in the Si 
slab reads 

Q = VgC^ {Vg - Vs)C, = {Vs - Vc)Csx = V,C s Z , (4) 

with the total capacitance C = [1/Ci + l/Csx + ^/Csz]^^ ■ 
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FIG. 2: (a) The charge density {p)y,z versus the x coordinate 
for model Si2- (b) The electrostatic potential {V}y,z versus 
the X coordinate for model S12. (c) The same as (b) for 
model L12. Curves at positive (thick) and negative (thin) gate 
voltages Vg of values 0.0, 4.35, 8.70, and 13.05 V are plotted 
as the dotted, solid, dashed, and dash-dotted respectively. 
The vertical dashed lines at a; = 2.86 and 8.77 A indicate the 
positions of the two surfaces between the device channel and 
the insulator slabs. The positions of atomic monolayers are 
denoted by the big (Si) and small (H) spheres. The average 
{. . . )y,z is done over y and z in the region < y < 5.43 A and 
8.15 <z < 13.58A. 

Here Vs is the surface potential and Vc the central poten- 
tial of the channel as further specified in Sec. III. B. 

A. Charge and Potential Distributions in Real 
Space 

The charge density {p)y,z distribution along the gate 
(x) direction is shown in Fig. [UJa) for model 812- Sim- 
ilar distribution is observed for other models and is not 
plotted here. Attracted by the charge in the metal gates, 
most of the net charges appear near the surface of the 
Si slab. As a result, as illustrated in Figs. [Sfb) and (c), 
the electric field concentrates mainly in the vacuum slab 
where {V)y^z varies rapidly and is weak inside the Si slab 
where {V)y^z becomes flat. The FET is in its off region 
at small Vg and {V)y^z shifts with Vg until Vg reaches the 
threshold voltage where {V)y^z saturates in the center of 
the Si slab. 

An interesting short gate effect is observed when com- 
paring the potential distribution for the short gate model 
S12 and that for the long gate model L12. For the 
short gate model S'12, as shown in Fig. [2Jb), {V)y^z ver- 
sus X in the vacuum regions (0 < x < 2.86A or 8.77A 
< a; < 11. 63 A) is nonlinear. This is a result of the fact 
that the length of the gate in model S'12 is comparable to 
the thickness of the vacuum slabs in model Si2- The elec- 
tric field in the vacuum slabs deviates greatly from along 
the X direction and the parallel-plate capacitor model 
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FIG. 3: The total net charge Q in (a), (c), and (e) and the 
surface potential Vs in (b), (d), and (f) versus the gate volt- 
age Vg are shown for models S\i (solid curve), 5i2 (dashed), 
521 (dotted), and I/12 (dash-dotted). Results for models of 
different insulator thickness are compared in (a) and (b), of 
different gate length in (c) and (d), and of different Si slab 
thickness in (e) and (f). The short bars beside the curves in 
(b), (d), and (f) indicate the Vs values at the threshold points 
where the Vs versus Vg curves change their slopes. 



fails. For longer gate model L12, the fringe effect on 
the channel region is reduced and the {V)y^z versus x 
curves in the vacuum slabs become straight as shown in 
Fig.HUc). The short gate effect described above is unique 
in nanoscale MISFET and does not appear in conven- 
tional MISFET since there the insulator slab's thickness 
is usually much smaller than its length. We will further 
discuss its effect on the FET's performance in Fig. [3jc) 
and (d). 



B. Total Transferred Charge and Surface Potential 

We use the surface just enclosing the H and Si atoms 
as indicated by the horizontal dashed lines in Fig. [IJa) or 
the vertical dashed lines at 2.86 A and 8.77 A for model 
S12 in Fig. [UJa) as the surfaces separating the channel 
and the insulator slab. The averaged value (^)y,z on 
these surfaces is used as the surface potential Vs. The 
central potential Vc is approximated as the {V)y^z value 
in the middle of the two gates, i.e. at a; = 5.82 A in 
Fig. Ha) for model S12. 

With Vs and Vc known at each Vg^ we can estimate 
the capacitance from Eq. For traditional MlSFETs, 
usually we have Csx < Ci in the off region and Csx > Ci 
in the on region. However, the Si slab here is extremely 
thin and Csx > Ci in both regions. In addition, the 
span of the device in the z direction is larger than that 
in the x direction. This fact leads to Csx > Csz in the 
off-region and Cs may increase with the Si slab thickness 
because Csz is proportional to the thickness while Csx is 



proportional to the inverse of it. In the on region, when 
a large amount of free charges are created on the surface, 
Vc is almost pinned due to the screening effect of these 
free charges. 

The total net charge Q versus the gate voltage Vg for 
model 5*12 is shown by the dashed curve in Fig. [3)^ a) . At 
Vg — the Fermi energy Ep is located in the energy 
gap of the Si slab. As Vg varies, the relative position 
of the Fermi level in the energy band shifts accordingly. 
Cs remains almost constant when the Fermi level is in 
the energy gap but increases quickly as the Fermi level 
enters the conduction or the valence band. This sud- 
den increase of Cs happens at the positive and negative 
threshold points, i.e., near Vg = 13 V and —7 V, respec- 
tively. Near these two points, due to the shift of the 
Cs/Ci ratio, the Vs versus Vg curve also changes its slope 
as shown by the dashed curve in Fig. [2Ib). 

To show the effect of the vacuum-slab thickness on 
the electrostatic characteristics, we use model with a 
thinner vacuum slab and redo the calculation. {p)y^z and 
{V)y^z have similar profiles between the gates as that for 
model 15*12 but with a larger amplitude at the same gate 
voltage. In the equivalent capacitive circuit, Ci increases 
with the shrinking of the vacuum slab, which results in 
the shift of Q and Vs as shown by the solid curve in 
Fig. (Sl^a) and (b). The Q versus Vg curve of thicker in- 
sulator slab model has a wider subthreshold region 
in the Vg axis. One interesting observation is that the 
threshold surface potentials, the surface potentials at the 
threshold points, for model 6*11 and S'12 are very close as 
indicated by the short bars beside the curves in Fig.[3l^b). 
This is because the surface potential reflects directly the 
Fermi level in the energy band of the Si slab. For a sys- 
tem of thicker insulator or smaller Ci it takes a bigger 
gate voltage difference to shift the Fermi energy from the 
valence band to the conduction band. 

In Fig. [2] we have shown that the electric field in the 
insulator slab can deviate from normal to the slab due to 
the fringe effect in short gate cases. This short gate ef- 
fect in nano FET reduces the per area capacitance of the 
insulator slab and results in a larger threshold gate volt- 
age than estimated from parallel capacitance approxima- 
tion. In Fig. ^c) and (d) we plot the Q versus Vg curves 
and the Vs versus Vg curves respectively for model S'12 
(dashed) and longer gate model L12 (dash-dotted). The 
smaller total capacitance C in model S'12, due to its half 
gate length plus the short gate effect in it, results in a 
less than half net charge Q in the system at any Vg and a 
much gentler slope of its Q versus Vg curve in Fig.[3l^c). In 
Fig. [3^d), we also observe a gentler slope of the Vs versus 
Vg curve for model S12 in the subthreshold region. Com- 
bining the relation Vs/Vg = Ci/{Ci + Cs) as expressed 
in Eq. ([U, this means a smaller ratio of CijCs and con- 
firms the conclusion drawn from Fig. [2] that the short 
gate effect mainly reduces the insulator capacitance. 

Up to now, we have seen that both making the insula- 
tor slab thinner and suppressing the short gate effect can 
enhance the insulator capacitance and reduce the thresh- 
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old gate voltage as illustrated in Fig. [3I^b) and (d). On 
the other hand, a close comparison of the curves shows 
that the two changes to the system have opposite effects 
on the threshold surface potential. The threshold sur- 
face potential has a tendency to increase in the former 
case (the solid bar on the positive Vg side in Fig. ^h) 
is slightly higher than the dashed one) while the thresh- 
old surface potential decreases in the latter case. A de- 
tailed study shows that the variation of the threshold 
surface potential is related to the potential distribution 
inside the Si slab. Given the same surface potential in 
the subthreshold region, the potential distribution along 
X changes little when varying the insulator slab thickness 
but a potential drop appears and broadens in the center 
of the {V) y,z versus x curves when shortening the gate. 

In Fig.[3l^e) and (f), we present the dependence of elec- 
trostatic characteristics on the Si slab thickness exempli- 
fied by model Sn and 5*21. Similar to the case of thicker 
insulator slab as shown in Fig. Efa) and (b) , the Q and Vs 
versus Vg curves for a system of thicker Si slab also have 
a wider subthreshold region in the Vg axis. However, this 
apparent common feature in the two cases has different 
origins and it is interesting to compare the details of their 
electrostatic characteristics. 

At first, as pointed out earlier in this section, because 
the Si slab inside the supercell has a dimension in the z di- 
rection bigger than that in the x direction, \Csx \ > \Csz\ 
and Cs increases with the Si slab thickness in the off- 
region. In this case, the Q versus Vg curve of thicker Si 
slab model S21 (dotted) in Fig. [3Ke) has a sharper slope 
in the off-region due to an increased total capacitance C, 
in contrast to the gentler slope of the curve for thicker in- 
sulator slab model (dashed) in Fig.[31[a). In addition, 
the bigger Cs in the thicker Si slab case, in stead of the 
smaller Ci as in the thicker insulator slab case, results 
in the gentler slope of the dashed Vg versus Vg curve in 
Fig.lSf). 

Secondly, it is well known [H, [13] that the energy 
gap of semiconductor materials in nanostructures is en- 
larged from their bulk value as a result of the quantum 
confinement effect. The change of the gap depends on 
the confining potential profile characterized by parame- 
ters such as the confining dimension and the potential 
height. This effect is automatically taken into account 
by our NEGF-DFT model. As will be illustrated later in 
Sec. III. C, the thicker Si slab in model S21 has a narrower 
and shifted energy gap compared to that in model Sn. 
As a result, the subthreshold region in the Vs axis, i.e. 
the difference between the surface potentials at the posi- 
tive and negative threshold points, of model 5*21 (3.79 V 
for the dotted curve in Fig. El^f)) is narrower than that 
of model Su (3.99 V for the solid curve). Furthermore, 
we also observe a shift of the subthreshold region in the 
Vs and Vg axes when the Si slab thickness varies. 

As discussed above and illustrated in Fig. [3Ke) and 
(f), the thicker Si slab in model 521, compared to that 
in model Su, narrows the subthreshold Vs region via the 
quantum confinement effect on the one hand and reduces 
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FIG. 4: (a) The transmission T and (b) the density of states 
(DOS) for electrons of momentum (0,0, fcz) are plotted as 
functions of the electron energy e in models Su (solid curve) 
and S21 (dashed). 



Vs/Vg ratio in the off-region due to the specific length- 
thickness ratio of the Si slabs in the models on the other 
hand. The observed wider subthreshold Vg region is a 
result of the competition between these two independent 
but opposite factors. 

In this section, we analyze the electrostatic character- 
istics of our model nano FET in a way widely employed in 
the analysis of conventional MISFET. [26j| Here we want 
to emphasize that our system (of size 2 nm) is much 
smaller than any of the conventional MISFETs even the 
smallest predicted by ITRS in its latest report. [l[ As 
a result, our model device shows characteristics between 
molecular electronics and conventional electronics. For 
example, we observe a large slope of the Q versus Vg 
curve in the off region in all models which is similar to 
the result of the Si4 cluster but different from that in 
conventional MISFET where a flat Q-Vg dependence is 
usually found in the off region. In addition, the threshold 
gate voltage in our systems (> 5 V) is much higher than 
the value usually observed in conventional MISFET (< 2 
V) . This is because the Si slab in our systems is very thin 
and a vacuum instead of Si02 slab is used as the insula- 
tor slab. The Cs/Ci ratio then becomes much larger in 
the off region. 



C. Transmission and Density of States 

To explore further the origin of the electrostatic prop- 
erties of nano FETs and the gate effect on transport, 
we calculate the transmission, which gives the conduc- 
tance in the linear region via Eq. ([3]) , and the density of 
states (DOS). In Fig. |D(a), the transmission of an elec- 
tron propagating along the z direction is calculated at 
= for quantum transport between the source and 
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FIG. 5: The transmission T versus energy e at |Vg| = 0.0 
(dotted), 4.35 (solid), 8.70 (dashed), and 13.05 V (dot- 
dashed) is illustrated for negative Vg (a) and positive Vg (b). 
Model Sii is used in the calculation. 

drain electrodes in model S\\ (solid) and S'21 (dashed). 
Because the electron sees a perfect crystal in the z di- 
rection, its transmission is 100% for each channel and 
the transmission versus energy curve tells the number 
of propagating channels for an electron of the energy. 
The one-dimensional DOS in fc^ direction for model 
(sohd) and S^x (dashed) is plotted in Fig. g^b). The 
band gap of model 5*1 1 is wider than that of model S'21 
due to the quantum-confinement effects in nano semicon- 
ductor structure [9J. In addition, most probably due to 
the change of the surface effect or the ratio between the 
numbers of the H and Si atoms in the system, a shift 
of the energy gap is observed when varying the Si slab 
thickness. 

The transmission versus electron energy at different 
gate voltages is plotted in Fig. [5] for model 6*11. When a 
negative gate voltage is applied, the energy band in the 
gate region shifts to higher energy. This shift forms an 
energy barrier for electrons and a well for holes. As a 
result, we observe a significant decrease of transmission 
for electrons near the bottom of the conduction band as 
shown in Fig. [5l^a). On the contrary, if a positive gate 
voltage is applied, an energy barrier is formed for the 
holes in the middle of the device and the propagation 
of holes near the top of the valence band becomes unfa- 
vored. The transmission profile changes quickly with the 
gate voltage in the subthreshold region and then becomes 
less sensitive when and Vs saturate at Vg higher than 
the threshold gate voltage. 

IV. SUMMARY 

In summary, we have proposed an atomistic multi- 
terminal model for electrostatic simulation of nano MIS- 
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FETs in the framework of the NEGF-DFT approach and 
studied the behaviour of devices with Si slabs of one and 
two unit cells in thickness, with variable insulator (vac- 
uum in this model) thickness, and with different gate 
length. The subthreshold region in terms of the surface 
potential reflects the band gap of the Si slabs and the 
surface potential at the threshold points remains little 
changed in devices of different insulator thickness. For 
a device with thicker Si slab, the subthreshold region in 
surface potential becomes narrower due to the quantum 
confinement effect while the subthreshold region in gate 
voltage becomes wider as a result of the competition be- 
tween the confinement effect and the geometry effect on 
capacitance. In nano FETs where the gate length is com- 
parable to the insulator thickness, the short gate effect 
results in a much wider subthreshold region than that 
observed in conventional MISFETs. This short gate ef- 
fect can be suppressed by using longer gates to reduce 
the fringe effect on the electric field in the insulator slab. 

In addition, we have calculated the band structure of Si 
slabs with different thickness. The density of states and 
the transmission probability of electrons in our nano FET 
system are then estimated to demonstrate the quantum 
confinement effect and the gating effect on transport in 
the linear region. The energy gap widens and shifts to 
higher energy when the Si slab varies from one to two 
unit-cells in thickness. The application of a gate voltage 
introduces energy barriers or wells to the carriers in the 
transport channel and modifies the performance the nano 
FETs. A negative gate voltage favors the hole transport 
between the electrodes while a positive one favors the 
electron transport. 

Limited to the computation capability, our model sys- 
tem is still not fully realistic and is smaller than even 
the latest feasible MISFETs [ij. Nevertheless, our model 
may provide qualitative information about realistic MIS- 
FETs and serve as a prototype framework for nano MIS- 
FET simulations. 
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